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We present the linear-stability analysis of synchronised 
states in coupled time-delay systems. There exists a synchro- 
nisation threshold, for which we derive upper bounds, which 
does not depend on the delay time. We prove that at least 
for scalar time-delay systems synchronisation is achieved by 
transmitting a single scalar signal, even if the synchronised 
solution is given by a high-dimensional chaotic state with a 
large number of positive Lyapunov-exponents. The analyt- 
ical results are compared with numerical simulations of two 
coupled Mackey-Glass equations. 
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The problem of synchronisation of dynamical systems 
is one of the classical fields in engineering science 0. 
Recently, renewed interest in this field was stimulated in 
connection with the synchronisation of chaotic motion. 
Especially, the potential applicability for communication 
has attracted much research in recent years 1^ . Yet, there 
are a lot of results available concerning the synchronisa- 
tion of low-dimensional chaotic systems, theoretical as 
well as experimental Contrary, the synchronisation 
of high-dimensional chaotic systems with possibly a large 
number of positive Lyapunov-exponents remains open. 
From the point of view of numerical simulations the syn- 
chronisation of specific high-dimensional chaotic systems 
has been achieved jsj, while, to our best knowledge, rig- 
orous results, e. g. concerning sufficient synchronisation 
conditions, are still lacking. While it has been proved 
recently by Stojanovski et al. Q that the synchronisa- 
tion of high-dimensional chaotic states can be in principle 
achieved with a single transmitted variable, the problem 
of finding the appropriate coupling of the two systems re- 
mains open. For that reason we address in this paper the 
question of the synchronisation of coupled identical time- 
delay systems. We focus on time-delay systems, since on 
the one hand it is well established that these systems are 
prominent examples of high-dimensional chaotic motion 
with a large number of positive Lyapunov-exponents , 
and one the other hand synchronisation of Mackey-Glass 
type electronic oscillators has been reported from the ex- 
perimental point of view |Q . 
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Let us consider a fairly general theoretical model and 
investigate the stability problem of a synchronised state. 
For that purpose consider two identical arbitrary scalar 
time-delay systems with a symmetric coupling 
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F{x,Xr) - K{x - y) 

F{y,yr) ~ K{y ~ x) 



(1) 



where we adopt the notation Xt ■= x(t — t) to indicate 
the time-delayed variables. We specialise from the begin- 
ning to the frequently analysed case that the coupling is 
bi-directional and acts additive to the single dynamical 
system. However, we stress that the subsequent consid- 
erations apply with minor modifications to much more 
general situations, e. g. to vector-type variables, to sys- 
tems with much more general delay terms, or to a non- 
additive coupling, as long as the coupling vanishes in the 
synchronised state x{t) = y{t). But we think, that the 
choice made in eq.(|l|) makes our arguments more trans- 
parent. 

Let z denote the synchronised solution, i.e. z — 
F{z,Zt). Considering deviations from that state ac- 
cording to a: = z -\- 5x, y = z + 5y and performing 
a linear stability analysis, we obtain for the deviation 
A :— Sy — Sx from the synchronised state the linear 
differential-difference equation 



A = a{t)A + /3{t)Ar 



(2) 
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Here, the time-dependent coefficients are given in terms 
of the synchronised solution as a{t) = diF{z, Zr) — 2K 
and f3{t) — d2F{z, Zr), where the symbol 9i/2 denotes 
the derivative with respect to the first /second argument. 
A superficial inspection of eq.(||) might suggest that the 
synchronised solution is stable if a{t) is "sufficiently neg- 
ative". In fact, we will make this statement rigorous in 
what follows. Suppose the coefficients are bounded in the 
sense that a{t) < —a < and \(3{t)\ < b holds for some 
fixed values a and b. Since the equation is linear, it is 
sufficient to analyse the solution with the special initial 
condition A(0) = 1, A{t) = 0,t < 0. The general case 
follows by a simple integration. There are different ways 
to estimate the stability of the trivial solution, A{t) = 0, 
of eq.(||). Here we use the fact that for scalar quantities 
a simple closed analytical formula for the solution can 
be written down. One just integrates the linear equation 
(^) in the time intervals [Nt, {N + 1)t] and considers the 
delay term as an inhomogeneous part. By this continu- 
ation method (cf. p. 45]) the full solution is obtained 
as 
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A(t) = e 



r ' a(t') dt' 
— P.Jo 



dt 
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+ f dti f dt2l3{h)l3{t2 - r)e4 "^"^ + 
+ dti dt2... dtN 
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, . , . , , ^ ^ f a(0)dB 

for TVr < t < (iV+ l)r . (3) 

Here, the domains of integration for the exponents are 
given by h := [0, ti - r] U [is - r, - 2r] U . . . U 

[tk — {k — 1)t, tfe — kr]). An upper bound for |A(i)| is 
obtained, if the maximal values a(t) = —a and (3{t) = b 
are inserted into eq.(^). But then, the expression reduces 
to a solution F of the differential-difference equation with 
constant coefficients 
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-aT + bV^ 



(4) 



Hence, a solution of eq.(Q) yields an upper bond for 
|A(i)|. But the last equation is easily solved by a Laplace 
transformation (cf. |^) or loosely speaking by an expo- 
nential ansatz T(t) = e*'*. Since the corresponding eigen- 
values obey s — —a + 6exp(— st), negative real parts, 
i. e. stability, occur if and only if a > 6. This inequal- 
ity yields an upper bond for the critical coupling 
strength beyond which synchronisation is achieved. If 
we take the definitions of a(t) and f3{t) into account it 
reads explicitely 



l/2[max9iF(z,z^) +max|a2F(z,z^)|] . (5) 



We note as a by-product that eq. (Q) may be viewed as a 
kind of Grownwall-like lemma |^ for the time-dependent 
equation (||). 

In what follows, we compare our analytical result to nu- 
merical simulations. We specialise to the Mackey-Glass 
system, i. e. 



F{x,Xt) 



1 



(6) 



In order to investigate the properties of the synchroni- 
sation mechanism by numerical methods, we chose the 
distance between trajectories as a suitable measure. For 
that reason the quantity 



Drit) 



t+T 



I xit') - y(t') I dt', 



(7) 



which of course depends on the range of averaging T and 
the point of reference t was analysed. 

We used a Runge-Kutta algorithm of fourth order with 
step size 0.1. The simulation, which have been performed 
for the parameter value a = 3, started with K = 0.35. 
A constant initial condition for x and y, which differs by 



an amount of 10~^ has been chosen. The system was 
allowed to relax for a time t = 80t. After that, the 
distance D was integrated on a trajectory of the length 
T = 80r. For the next value of coupling strength K we 
again distorted the last state of the trajectory by adding 
an amount of 10^"^ to the y-coordinates and used it as 
initial condition. We performed the computation for in- 
creasing as well as decreasing coupling constant. Fig.|l| 
summarises our findings. 

FIG. 1. Distance D for two coupled Mackey-Glass systems 
for r = 10 (dashed line) and r = 100 (solid line). 



FIG. 2. Time series for r = 10 and different values of the 
coupling (a) K = 0.2595, (b) K = 0.250, (c) K = 0.240. 

For r = 10 we observe distinct jumps, indicating that 
the system switches between coexisting periodic states 
with a pronounced hysteresis. For t = 100, the overall 
behaviour of the system appears to be quite similar as 
for T = 10, except that no switching and no hysteresis is 
observed. Within the resolution of the graphics the same 
behaviour has already been observed for a smaller value 
T — bQ. From the numerical simulations the synchro- 
nisation threshold is estimated as Kc{t — 10) w 0.24, 
Kc{t — 100) « 0.28. If we evaluate our analytical 
estimate eqs.(^ and using upper bounds for the 
derivatives we obtain values which differ by an order 
of magnitude but are independent of the delay time r, 
K+ = (81a/40 - l)/2 = 2.5375. Since we have applied 
a rather graceful rigorous estimate such a discrepancy is 
far from being astonishing. 

In order to understand the dynamics in the vicinity of 
the synchronisation threshold, time traces of the differ- 
ence x{t) — y{t) have been computed (cf. fig.|l). Slightly 
below the synchronisation threshold Kc we observe an 
intermittent behaviour very similar to on-off intermit- 
tency |l0|] . Additionally, we investigated the distribution 
of laminar phases and turbulent phases under variation 
of the coupling strength K and the delay time r, which 
we present in fig. ||. To this end, the distance Dr(t) 
of the two systems in the phase space has been com- 
puted on an intervall of length rlO^. Then, the length 
of the laminar phases {Dr < 0.10) and turbulent phases 
{Dt > 0.10) were recorded. We observe a power-law 
scaling of the distribution Pi of the laminar phases over 
a wide range, P; oc t~'^'^'^\ where the exponent ai de- 
pends slightly on the delay time. In the low-dimensional 
chaotic case, for r — 10.0, we observe a(10) = 1.50 in 
agreement with the value predicted by the scenario of on- 
off intermittency. In the high-dimensional chaotic cases, 
for increasing delay time, we observe a decreasing expo- 
nent: a(30.0) = 1.41,a(50.0) = 1.38, a(lOO.O) = 1.27, 
indication that there might be deviations from on-off in- 
termittency. The distribution of the turbulent phases 
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follow a Pf oc i "'("^^-scaling for high enough r. 

FIG. 3. (a) Distribution of laminar phases, (b) distribution 
of turbulent phases for r — 100. (solid line: K = Kc ~ 0.28; 
circles: K = 0.26, dotted line: K = 0.24, dashed line: 
shifted power-law fit for K = Kc ~ 0.28 with q; — 1.27, 
and Qt = 2.19 

Although these findings are at the first sight not sur- 
prising, a detailed analysis is required, with special focus 
on the high dimensionality of the dynamics. Details will 
be reported elsewhere. 

We conclude with a remark, how our results depend 
on noise or other imperfections which are present in re- 
alistic systems. In fact, in order to apply a concept like 
synchronisation such perturbations have to be small and 
we may assume a general linear dependence. Formally 
such contributions are introduced into eq.(^ by adding 
the two terms G{x,Xt)£, and G{y,yr)r], where ^ and r] 
denote for example realisations of a noise. Considering 
the perturbations of the same order of magnitude like the 
deviation from the unperturbed synchronised state and 
proceeding as above we finally end up with 

A = a{t)A + (3{t)^r + G{z,Zr){v~0 , (8) 

which differs from eq.(|^) just by an inhomogeneous con- 
tribution. The theory of linear difference-differential 
equations tells us Q that eq. (^) inherits its stability prop- 
erties from the corresponding homogeneous system (|^) 
except that the perturbations cause fluctuations around 
the unperturbed synchronised state. Whenever the per- 
turbations are so large, that contributions beyond the 
linear order have to be taken into account, one has to 
resort to different methods. One of these cases, which 
are also relevant from the experimental point of view, 
is given by the synchronisation of nearly-identical time- 
delay systems. Since, in this case, no strict synchronised 
solution X = y exists, one has to rely on more general 
concepts, such as the generalised synchronisation [pi]] . 

In summary, we emphasise that an analytical upper 
bond for the solution of eq.(||) is obtained if one replaces 
the time dependency of the coefficients by their extreme 
values. One might get better estimates in special cases. 
In particular one might argue, that eq.(Q) already deter- 
mines the stability, if the time averages of the coefficients 
are inserted. This statement is in fact true if either the 
coefficients are periodic functions of time with the delay 
T being an integer multiple of the period or, if the co- 
efficients are almost constant (cf. p. 277]). Whether 
the general case can be treated by this refined estimate 
remains open. Nevertheless we have shown, that for suf- 
ficiently large coupling constant K the synchronised so- 
lution of eq.(|l|) becomes stable, whenever \5i/2F{z^Zr)\ 
are uniformly bounded. In particular the critical cou- 
pling strength remains bounded even in the limit of large 
delay times, i. e. it does not increase with the dimen- 



sion of the attractor. In fact, our numerical simulations 
indicate only a weak dependence of the actual critical 
coupling strength on the delay time. Last but not least 
our approach clearly demonstrates that the success of the 
synchronisation is independent of the number of positive 
Lyapunov-exponents, even if our coupling uses one scalar 
variable only, illustrating the results of Stojanovski et al. 
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